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The effective interactions of dark matter with photons are fairly restricted. Yet both direct 
detection as well as monochromatic j ray signatures depend sensitively on the presence of such 
interactions. For a Dirac fermion, electromagnetic dipoles are possible, but are very constrained. 
For Majorana fermions, no such terms are allowed. We consider signals of an effective theory with 
a Majorana dark matter particle and its couplings to photons. In the presence of a nearby excited 
state, there is the possibility of a magnetic dipole transition (Magnetic inelastic Dark Matter or 
MiDM), which yields both direct and indirect detection signals, and, intriguingly, yields essentially 
the same size over a wide range of dipole strengths. Absent an excited state, the leading interaction 
of WIMPs is similar to the Rayleigh scattering of low energy photons from neutral atoms, which 
may be captured by an effective operator of dimension 7 of the form XX.Fp.vF^ . While it can be 
thought of as a phase of the Magnetic inelastic Dark Matter scenario where the excited state is 
much heavier than the ground state, it can arise from other theories as well. We study the resulting 
phenomenology of this scenario: gamma ray lines from the annihilation of WIMPs; nuclear recoils 
in direct detection; and direct production of the WIMP pair in high-energy colliders. Considering 
recent evidence in particular for a 130 GeV line from the galactic center, we discuss the detection 
prospects at upcoming experiments. 



I. INTRODUCTION 

Weakly interacting massive particles (WIMPs) have 
long been studied as potential candidates for the cold 
dark matter observed in the Universe. The most well- 
motivated and deservedly most well-studied WIMPs are 
those that emerge in extensions of the Standard Model 
associated with the seemingly unrelated problems of the 
electroweak scale, such as supersymmetric extensions. 
An orthogonal line of inquiry is motivated by the de- 
ceptively elementary question of "how dark is Dark Mat- 
ter?" Namely, what are the strongest constraints on the 
interaction of dark matter with the electromagnetic field? 
Numerous studies already exist and in particular the idea 
of electric and magnetic dipole interactions have recently 
attracted considerable attention [TMll|. In these models, 
single photon exchange provides a possible direct detec- 
tion signal, while annihilation into two photons might 
provide an indirect detection signal (see e.g., [T2"]). 

However, if dark matter is a Majorana fermion, then 
these single-photon couplings through electromagnetic 
dipoles do not exist - the dipole operator vanishes identi- 
cally for Majorana fermions. For Dirac fermions, it is nat- 
urally off-diagonal [131 E]. For pseudo-Dirac fermions, 
in which case the ground state is the dark matter candi- 
date, the dipole interaction mediates transitions between 
this ground state x and an excited state x* ■ The authors 
of Ref. [5] exploited this possibility to build a model, 
dubbed Magnetic inelastic Dark Matter (MiDM), to ex- 
plain the DAMA results through dipole-dipole dominated 
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scattering. The interaction Lagrangian of MiDM is 

where /i x is the dipole strength, B^ u is the hypercharge 
field-strength tensor, and <7 P „ = £[7^, 7„]/2. This cou- 
pling contains within it the interaction with the electro- 
magnetic field. We study the signatures of this model in 
the first part of this paper. 

In the limit that we take the excited state heavy, we 
are left with a Majorana fermion, and we can again ask 
the question "how dark is Dark Matter?" . Starting with 
the MiDM Lagrangian above if the excited state, x* > is 
much heavier than the energy available then it can be 
integrated out to yield the interactions 

2 

C = (XXB^B^ + ixi.xB^B^) , (2) 

X* 

where W v = \e^ af3 B a0 and e'" /a ' ? is the Levi-Civita 
symbol. Motivated by this form, in the second part of this 
paper we will concentrate on the slightly more general 
case for the interaction of DM with the electroweak field 
strengths 

C =lhr \xx (cos OB^B^ + sin O^W^W^) (3) 
+ i X75X (cos x B^B^ + sm8 x TvW^W^ } . 

Here quantifies the relative coupling to the field 
strength of hypercharge in comparison to that of SU W (2) 
and A R is some high scale related to the cut-off scale of 
the theory. We will discuss UV realizations in a later 
section, but simple scenarios can arise either as a limit of 
MiDM, or for instance integrating out a dilaton (or axi- 
dilaton). The interactions of Eq. ^ are akin to the fa- 
miliar interactions of photons with neutral atoms at long 
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wavelengths that lead to Rayleigh scattering. Hence we 
dub this scenario Rayleigh Dark Matter (RayDM) . This 
could be the entirety of the DM interaction with the stan- 
dard model, but it also serves as a reasonable form of the 
effective operators responsible for 7 lines in many mod- 
els (even when they freeze out dominantly through other 
channels). The special form of this interaction, which 
necessitates at least two force mediators, requires a re- 
consideration of the basic processes by which we hope 
to detect dark matter and this constitutes a part of the 
current work. 

In this paper we set to explore these different possibil- 
ities for the interaction of Majorana WIMPS with light. 
The paper is organized as follows: In section [TT| we dis- 
cuss in detail the MiDM scenario including its signatures 
in gamma rays as well as the prospects for seeing it in 
direct detection experiments; In sectio n |III| we explore 
the phenomenology of RayDM; Section |IV| is devoted to 
the prospects of collider searches for both MiDM as well 
as RayDM; Finally, the main findings of this work are 
summarized in the conclusions, section [VJ We caution 
the reader that the clear separation between MiDM and 
RayDM is not always appropriate. As we shall discuss 
and emphasize below, there are certain aspects of the 
phenomenology where the two scenarios and the opera- 
tors involved cannot be logically separated. 



abundance is determined by the annihilation into photons 
and vector-bosons, which typically requires larger dipole 
strength. The prospects for direct detection in this case 
are fairly gloomy, but collider constraints provide strong 
and interesting bounds on this possibility as we describe 
in section IIVI 

As a consequence, in this section we focus on the case 
when the mass splitting is small. In particular, when the 
mass splitting to the excited state vanishes or is smaller 
than the kinetic energy of the WIMP in the halo AM = 
m . — m ~ 100 keV the cross-section of scattering on 
the nucleus is much larger and the corresponding rates 
in direct detection experiments are phenomenologically 
interesting. 

Before moving onto the details, it is important to 
make the phenomenological point: for thermally produced 
MiDM, the direct and gamma-ray line indirect signatures 
are roughly independent of the size of the dipol¥\ That 
is to say, even if MiDM is a sub-dominant component of 
the dark matter, the amplitude of these signals would be 
unchanged. 

This arises quite simply. The annihilation rate is 
governed by the annihilation of WIMPs into fermion 
pairs |15j which scales as /i 2 . This implies that the num- 



ber density of WIMPs scales as n x 



II. MiDM 



w /^thermal 

PMiDM — PO x n : 

Px 



(5) 



In this section we concentrate on the MiDM scenario, 
but consider a slightly more general form of the magnetic 
dipole interactions 
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where and are respectively the field strength 
of the photon and Z boson, u and /i z are the corre- 
sponding dipole strength, and x an d X* are two Weyl 
fermions. If the interaction above the electroweak scale 
is entirely with the field strength of hypercharge then 
Mz/A 1 -,. = — tan# w . In what follows we explore the more 
general possibility since one can entertain additional op- 
erators with the field strength of SU W (2) that would re- 
sult in a different ratio. An example of such an opera- 
tor is the dimension 7 operator (-^h) ^cr^xTr h^W^h. 
However, considering that such operators are generically 
further suppressed we expect the deviation away from 
the relation fJ, z /p> = — t&n9 w to be small. 

The phenomenology of this theory depends crucially 
on the mass splitting between the two states, x an d X* ■ 
When the mass splitting is large m x , — m x > m x the 
phenomenology is similar to that of RayDM, which is 
described in the next section. More generally, even for 
a smaller splitting such that m » > m x /20 the heavier 
state x* is n °t present during the early universe freeze- 
out of the lighter state %. In this case, unless other 
channels are available through new interactions, the relic 



where /^thermal is the dipole needed to achieve the ap- 
propriate relic abundance that arises from a cross sec- 



tion of 6 x 10 26 cm 3 s 



For a WIMP of mass m 



130 GeV we find thermal « 1.2 x W- 3 n N (2.0 x 10- 3 Mjv ) 
for a dipole ratio of /i z /u — — tan 6 W (/i z /a = 0) . 
Here p, N — 0.16 GeV^ 1 is the nuclear magneton. This 
dipole strength corresponds to an annihilation rate into 
gamma rays of cr(xx ~^ ll) = 2.5 x 10~ 29 cm 3 /s 
(1.6 x 10~ 28 cm 3 /s ). 

The collision rate in direct detection experiments scales 
as n x <7 X N, thus 



R DD cx n x p? x 



PO ^thermal v ,,2 

X [j, x 



P0 



/4 



Mthermal i (6) 



where po « 0.4 GeV/cm 3 is the local density of a WIMP 
comprising all of the dark matter. Such a scaling phe- 
nomenon is well known in many WIMP models, that 
when s-channel annihilation diagrams are directly linked 
to t-channel scattering, the lower relic abundance is com- 
pensated by the higher scattering cross section (see, e.g., 



1 Since the breaking of hypercharge makes the phenomenology de- 
pends in principle on both fj, as well as fi z , this statement as- 
sumes the absence of any unexpectedly large difference in scale 
between /t z and u . 
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A possibly more remarkable scaling is associated with 
the 77 and 7Z signatures. The cross section for these 
processes is proportional to /i 4 ,. Thus, the indirect 7-ray 
rate scales as 



Po ^thermal v ,,4 
9 A X Mv 



^L 4 

_9 /^thermal" 

X 



(7) 



Again, if the WIMP is thermal, the 77 signal is inde- 
pendent of the size of the dipole, even if the fraction of 
the dark matter is much smaller. There are important 
caveats to this, as we shall discuss below, but they do not 
change the fact that even for very large dipoles (yielding 
under-abundant dark matter) the signals are at the same 
level as that of a thermal dominant WIMF0 

However, such a scenario is excluded unless 8 — ra x * — 
m x > l/2m x v 2 and the direct detection scattering is ei- 
ther inelastic (if the excited state is accessible) or not 
present (if it is not). Intriguingly, the size of the signals 
would be appropriate for DAMA [T7| (via the MiDM sce- 
nario [5]) and approximately for the recently reported 130 
GeV signal, given the astrophysical uncertainties [TBH2T] . 



A. Annihilation Rate and Relic Abundance 

To understand these points in detail, let us consider 
the precise values realized. When the excited state x* 
is present, the annihilation into Standard Model fermion 
pairs through 7/Z dominates. The annihilation cross- 
section to leading order in velocity is given by 



eKxx* -> ff) v 



2 2 



(8) 



+ (*? 



Here qf is the fermion's electric charge, Vf (a/) is the 
ratio of its vector (axial-vector) coupling to the Z boson 
to its electromagnetic coupling, and £(s) = s/(s — m|). 
In the above expression we took = m x since it is 
only when the mass splitting is not too large that the 
heavier state is relevant. The annihilation rates of the 



lighter state into vector-bosons are 
111 ' / 2ir \ 2 



ct(xx* -> 7Z)« = 2 



<?(xx* 



ZZ)v = 




2m m 



I 1 + 



In the case of large mass splittings, we are effectively left 
with a single species at freeze-out and its annihilation 
into vector-bosons must therefore be sufficiently large so 
as to yield er tot u ss 3 x 10 _26 cm 3 /s. This implies a rather 
large dipole strength which is in tension with collider 
searches for mono-photons discussed in section |IV| 

In the case of intermediate and small mass splittings, 
the cosmological history is slightly different and the nec- 
essary total annihilation rate is consequently altered. 
Since the heavier state decays to the lighter state through 
the dipole transition only after freeze-out, we effectively 
have two species during freeze-out, each of which can 
only annihilate on the other. The relic abundance neces- 
sary at freeze-out is therefore only half its usual valu^J 
This implies that the total annihilation rate is larger 
CTtotW ~ 6 x lCP 26 cm 3 /s. Interestingly, since this requires 
a dipole strength larger by a factor of \/2, this leads to an 
increase of the annihilation rate into 77 and 7Z of factor 
by 4. 

In two recent papers, Refs. jTSJ [TS] reported on a 
tentative gamma ray line in the Fermi/LAT data at 
Ej = 130 GeV, which has recently been confirmed by 
Tempel et al. and Su and Finkbeiner [21] • This result 
can be accommodated within the MiDM scenario with 
a WIMP mass m x = 130 GeV leading to the gamma 
ray line at E 1 = 130 GeV through XX ~ * 77- It will 
generically also result in an additional line of compara- 
ble strength at E 1 = m x — m|/4m x =114 GeV from 
XX 7Z, which is consistent with the data (see also the 
discussion in ref. [22] J] In Fig. [l] we plot the annihilation 



2 We do not consider the implications of a CP violating inelas- 
tic electric dipole moment here. Due to a velocity-unsuppressed 
dipole-charge scattering (see e.g., [2]), the constraints in 1151 con- 
strain the dipole fi x < 10 _8 /^ N . At these levels the indirect sig- 
nals would be negligible, unless the excited state is inaccessible. 

s Once fj, x is large enough this scaling ceases because the annihi- 
lation to gauge bosons dominates. For 130 GeV this occurs at 
fi x > .05/^jv- At this size, we shall see that collider constraints 
would exclude the scenario already. 



4 This is true only as long as the heavier state's lifetime is suffi- 
ciently short so that present day DM is entirely composed of the 
lighter state. For m » — m v > 100 keV the lifetime is shorter 
than about a microsecond. 

J Alternatively, m x = 144 GeV may give rise to the line at 
E~ = 130 GeV through the annihilation \X ~ * 7^- In this case 
the annihilation to two photons would have to be somewhat sup- 
pressed since no feature is observed at Ej = 144 GeV. Such a 
suppression is more natural in the RayDM scenario discussed in 
the next section. 
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FIG. 1. The total annihilation rate for m„ = m « = 130 GeV 
to fermion pairs, a{xx ~^ 77) m the MiDM scenario as a 
function of the relative strength of the photon dipole to Z 
dipole. The dipole itself is normalized to yield the annihilation 
rate into 77 (blue band) and 7Z (green band). The bands 
span the range shown in units of 10 _27 cm 3 /s. The red vertical 
line shows the relevant ratio of dipoles in the case of MiDM 
interaction with hypercharge only. 

rate into fermions as a function of the dipole ratio /x z //u, 
for a WIMP mass of m = 130 GeV in the case of small 
splitting m x w m x , . 

Surprisingly, when the dipole strength is normalized 
to yield an annihilation rate into 77 (7Z) in the range 
recently reported in refs. QUOIT], 0.3-1.3 x 10" 27 cm 3 /s, 
the annihilation rate at freeze-out into // pairs near 
the expected ratio of /x z / fi = — tan 9 W is only a fac- 
tor of 3 — 7 (5 — 10) larger than the needed value of 
6 x 10~ 26 cm 3 /s. As we shall see in the next subsec- 
tion this also yields a rate in direct detection experi- 
ments which can be probed with existing experiments 
and is surprisingly close to that reported by the DAMA 
collaboration for m x — m . « 100 keV. Evidently, this 
surprising concordance is numerically not perfect, but is 
sufficiently interesting given the large astrophysical un- 
certainties. 



B. Direct Detection of MiDM: Constraining 
Cosmic Ray Gammas Underground 

A WIMP with a magnetic dipole can scatter against a 
nucleus' own magnetic dipole as well as its charge. In the 
case of elastic scattering the current limits from direct de- 
tection experiments on a WIMP of mass m x « 130 GeV 
and local mass density of 0.3 GeV/cm 3 are at the level of 
P"X £6x 10 _5 /i N (see e.g. [15]). This excludes annihila- 
tion rates into di-photons at the level of a(xx ~^ ll)v < 
10~ 34 cm 3 /s, far below anything we can hope to mea- 
sure anytime soon. Thus, current direct detection limits 
robustly exclude the possibility of observing gamma ray 
lines from magnetic dipoles annihilations independently 
of the dipole strength. 



A more promising possibility is offered by MiDM [5] 
where the WIMP couples via a magnetic dipole to an 
excited state x* ■ If the mass splitting is of order the ki- 
netic energy of the WIMP in the halo AM = m » — m x ss 
100 keV then it may undergo inelastic scattering against 
the nucleus, but the corresponding rates in direct detec- 
tion experiments are much reduced compared to the elas- 
tic scattering discussed in the last paragraph. Thus, the 
dipole strength can be larger and the annihilation rates 
of WIMPs into gamma rays considerably enhanced. This 
scenario then offers an interesting connection between 
observations in gamma ray lines and direct detection ef- 
forts that can be probed with current experiments. In 
particular, as discussed in ref. [9] it may explain the sig- 
nal claimed by the DAMA collaboration^] [T7] for WIMP 
masses of m x ~ 100 GeV and magnetic dipole strength 
in the range fi x = 10~ 2 — 10~ 3 /x N . Interestingly, for m x = 
130 GeV and fi x = 3 x 10~ 3 /j, N the annihilation rate of 
WIMPs into gamma rays is av(xx ~^ 77) ~ 10 - cm 3 /s 
which can accommodate the excess recently reported in 
ref. 0i. 

C. Variations on a theme: Model Dependences in 
Indirect Signals 

As we laid out in eq. [7] the 7-ray line signals should 
be independent of the size of the dipole. The question 
arises as to how robust the size of the gamma ray signal 
will be to changes in the underlying model. Indeed, we 
have already seen examples of this: an MiDM model that 
only coupled to 7 would be a thermal relic with [i x = 
^thermal ~ 2 x lO -3 /^ , which leads to a 77 signal for a 
signal size of (cry) ~ 1.6 x 10 _28 cm 3 s _1 (normalized to 
the case with Q x = Qdm)- I n contrast, in the presence 
of a dipole with the hypercharge gauge boson, the needed 
dipole for a thermal relic is roughly v3 smaller, leading 
to an effective signal size of (av) ~ 2.5 x 10 _29 cm 3 s _1 . 
While the pure 7 dipole is close enough to explain the 
signal at 130 GeV, a dipole of hypercharge seems too 
small except in very cuspy halos. Thus, we should inquire 
whether there are any effects that modify this. As it 
turns out, there are at least two simple elements without 
enlarging the effective theory that can affect this. 

The first effect is the presence of form factors. Since 
the appropriate dipole scale is fi x > e/ TeV, it should 
be resolved near the WIMP mass scale. The anni- 
hilation into 77 samples a form factor with spacelike 
q 2 = —rn 2 ,, while the annihilation into charged pairs sam- 
ples a form factor with timelike q 2 = Am 2 . (The 7Z 



6 Recently, the KIMS collaboration has claimed to exclude the pos- 
sibility by O(l) at 90% confidence |23l . However, this was a par- 
ticular range of energies arising for a specific choice of quenching 
factors, both on Nal and Csl. When combined with the absence 
of a thorough discussion of energy resolutions (given that sig- 
nal may leak into surrounding bins), the MiDM scenario appears 
intact, although likely requires O(l) modulation fraction. 
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FIG. 2. (left) Loop diagram contributing to dipole operator for MiDM model, (right) Comparable diagram contributing the 
Rayleigh operator. 



signal samples form factors with q 2 = ±(m 2 — mf/2).) 
If the form factor is being resolved at these momen- 
tum transfers, then treating it as a contact interaction 
is clearly wrong. Since q 2 for the s-channel diagram 
is four times larger than the t-channel, it is reason- 
able that it could pick up a larger suppression factor. 
In this case the dipole would have to be increased by 

Mx = Mthormal/ F2 ( 4m x)- The reSultin g si g nal in 77 

would be increased by (F 2 (—m 2 ,) / F 2 (4m 2 ,)) 2 . We note 
that in this case, if the scale A of new physics is not 
much higher than m x , the Rayleigh operator we discuss 
in section [III] may be generated as well, which could con- 
structively interfere, enhancing the 77 rate. 

It is worth dwelling on this last point for a mo- 
ment. Generally, we should in fact expect the presence 
of the Rayleigh operator whenever the dipole operator is 
present. If the dipole is generated by a loop process as 
shown in Fig. [2^, then assuming a coupling \\XX , the 
natural size for the dipole operator is 



A 2 



1 



9Y 



16tt 2 M x 



(12) 



By attaching a second external photon, the Rayleigh op- 
erator is also generated through the diagram shown in 
Fig. [2b with a natural size 



that the Rayleigh operator could contribute, it should be 
a likely expectation. 

A second possibility is the presence of CP violation. 
While both CP-conserving magnetic and CP-violating 
electric dipole moments produce s-wave 77 signals, only 
the s-channel diagram via a magnetic dipole yields an 
s-wave annihilation to charged fermions [15j . Thus, 
with EDMs one can increase the present day 77 sig- 
nal while producing only a p-wave suppressed annihila- 
tion at freezeout into //. With only EDMs, assuming 
a freezeout at T m x /20, the annihilation into gauge 
bosons dominates, and one has a signal cross section of 

6 x 1CP 26 cm 3 s _1 (as it is a Dirac fermion at freeze- 
out), exceeding the Fermi limits. Thus, while pure EDM 
is excluded, a combination of EDM and MDM could pro- 
duce the signal consistent with constraints. However, as 
we have stated, in the presence of such a large EDM, the 
direct detection limits would have excluded it unless the 
excited state is completely inaccessible. So while the first 
possibility still offer the prospect of discovery at upcom- 
ing direct detection experiments, this second case seems 
unlikely to be found underground. 



A 2 



1 



9y 16tt 2 M| 



(13) 



III. RayDM 



Thus, annihilation to 77 through the RayDM process is a 
one-loop process, while annihilation through the MiDM 2 
process is effectively two-loop. Thus RayDM annihilation 
is relatively enhanced over MiDM 2 by a factor of 



A 2 



16tt 2 



" 2 m£ 
Ml 



(14) 



Even for m x ~ Mx, if the theory is at all perturba- 
tive, the Rayleigh contribution should dominate. If the 
Rayleigh contribution to the amplitude is even a few 
times larger than the MiDM 2 contribution, the size of 
the signal should be easily large enough to explain the 
130 GeV signal. Thus, not only is it reasonable to believe 



A. Searching for RayDM in Gamma Rays 



The non-relativistic annihilation cross-section of 
RayDM into the different electroweak vector-bosons is 
sensitive only to the axial X75X components to leading 
order in the velocity expansion. The differential cross- 
section is given by 



ct(xx -> VV)v 



9 4 

47T A 6 



a; 



(15) 
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with the kinematic functions K. vv and couplings g vv de- 
fined as 



/Cey/v ' 1 - 



/C 7 z = 2 1 - 



K m = 1 



4m 2 



#77 = C x C w 2 + S x S w 2 > ( 16 ) 
5 7 Z = ( S x - C x) S w C w:( 17 ) 



3/2 



i 5zz = C x S w +Sx C w . (18) 



3/2 

^"WW = ^ ( 1 J , <?ww = S X ' 




(19) 



Here c = cos 9 and s = sin 9 and c,„ and s,„ are 

X XX x w w 

similarly defined with respect to the Weinberg angles. 
When m is not too much smaller than A R we expect 
some form-factor suppression to soften the behavior of 
this cross-section. 

In Fig. [3] we plot the annihilation cross-section of 
XX — > 77 as a function of the WIMP mass for several val- 
ues of cos 9 x , the relative coupling to the field-strengths 
in Eq. (j3j). Requiring the right relic abundance, which for 
a Majorana fermion is obtained when the total annihila- 
tion cross-section at freeze-out is 3 x 10 _26 cm 3 /s, we can 
normalize the Rayleigh scale A R . For m = 130 GeV 
this results in A R = 440 GeV (A R = 490 GeV) in the 
case of cos6> x = 1 (cos# x = 0). For a Dirac fermion 
the necessary annihilation cross-section is 6x 10 _26 cm 3 /s 
and the Rayleigh scale is correspondingly a factor of 2 1 / 6 
smaller. The resulting gamma rays are monochromatic 
with E~f — m x . To qualitatively understand these re- 
sults, we consider the limit where the WIMP mass is 
much heavier than the vector-bosons's. Then the expres- 
sion for the total cross-section is particularly simple and 
by equating it to the required cross-section from relic 
abundance we can solve for A R in terms of the WIMP 
mass and the angle 9 , 



X)°"(xX-> VV')v = 3 x 10~ 26 cm 3 /s 



(20) 



vv 



A = 600 GeV 



/ m \ 2 / 3 , 

(2^) (*-«»2* x ). 



With this value of the Rayleigh scale the annihilation 
rates into 77 and 7Z are 



^(XX -> 77)^ _ ( c w 2 c x + s w \) 2 ( K\ 



3 x 10- 26 cm 3 s- 1 2 -cos 29 



A r 



\v(XX ->■ 7^)f 
3 x 10- 26 cm 3 s- 1 



c 2 s 2 



(c -s ) 2 /A thN 

\ X X / J R 



2 - cos 29 



,(21) 



.(22) 



Here A' h is the value of the Rayleigh scale that leads to 
the correct relic abundance. As can be expected when the 
Rayleigh operator is mostly associated with hypercharge 
the total cross-section is very close to the cross-section 
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FIG. 3. The annihilation rate of WIMPs to di-photons, 
o(xX ~* 77), as a function of the angle 9 X for different choices 
of the WIMP mass. For each mass choice the Rayleigh scale 
A H is chosen so that the total annihilation cross-section yields 
the correct relic abundance. Shown are m x = 100 GeV (solid- 
black), m = 130 GeV (dashed-red). The dotted-blue curve 



depicts the asymptotic formula Eq. (21 1 which is independent 
of mass. 



for annihilation into a photon pair <7 t ot ~ c(xx — > 77), 
which would result in too large a signaF] When the 
Rayleigh operator is mostly associated with the non- 
abelian SU W (2) part, the annihilation into photons is 
suppressed compared with the total cross-section due 
to the Weinberg angle and the XX W + W~ channel. 
More quantitatively, in the pure SU W (2) case, <J gg /<Jtot ~ 
1/28 and 1/2 x a-yz/^tot ~ 1/12. The photon-to-hadron 
ratio (<r 77 + l/2a 7 z) /(Ptot — C77 — l/2c 7 z) ~ 1/7.6. In 
the case of hypercharge RayDM the equivalent numbers 
are a gg ja tot w 1/1.4, 1/2 x a lZ /^tot ~ 1/7. The photon- 
to-hadron ratio is w 5, so no significant hadronic emission 
is present in this case. 

An additional process contributing to monochromatic 
gamma ray signal is of course XX ~^ 7^ with a lower 
energy of = m x — to|/4to x . In Fig. [4] we plot the 
annihilation rate associated with this channel as well as 
its ratio to the di-photon rate. Depending on the DM 
halo profile and the angle cos 9 , both the 77 and 7Z 
rates are in the interesting range reported in ref. |19j . 
3 x 10~ 28 - 2 x 10 _27 cm 3 /s. This points to a fairly low 
Rayleigh scale of A R w 500 GeV. 

One might worry about the validity of this picture 
given that the Rayleigh scale is rather low. We come back 
to this point in section |TV| where this issue is particularly 
important, however, for the purpose of non-relativistic 
annihilation it is only necessary for the Rayleigh scale 
to be somewhat larger than the WIMP mass. Never- 
theless, since the WIMP mass is not much lower than 
the Rayleigh scale, it may be appropriate to include a 



7 Although, as we discuss in the conclusions, a subdominant DM 
component would plausibly give the right signal. 
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FIG. 4. The top pane is similar to Fig. [3] above, but for half 
the annihilation rate into a photon and a Z boson, 2 a (xX ~~ ^ 
7Z). In the bottom pane we plot the ratio of the annihilation 
rate into 7Z to that into 77 for m x = 130 GeV as a function 
of the angle 8 . 



form-factor. Thinking about RayDM as MiDM 2 allows 
to resolve the 4-point interaction with the exchange of 
the excited state x* ■ Consulting the corresponding an- 
nihilation rates in MiDM, eqs. (9pT ) we see for example 
that the process XX ~~ ^ 11 IS diminished by a factor of 



(1 



B. Direct Detection of RayDM 



The scattering of RayDM against matter is compli- 
cated by the fact that at least two force mediators have 
to be exchanged. The exchange of two photons leads to 
the least amount of suppression and so we concentrate 
on this case and consider the Lagrangian 
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c 



4A 3 „ 



XX F^F^ 



(23) 



Two distinct processes are possible: %A — > xA elastic 
scattering through the loop shown in Fig. [5} \N — > xA + 
7 tree-level scattering. The latter channel is extremely 
suppressed due to phase-space. The total cross-section 
for the elastic channel was previously calculated in the 
thorough work of ref. [5] in the approximation that the 
nucleus is much heavier than the WIMP by considering 



k' 



l + q/2 




P = v + l + Q /2 



FIG. 5. Elastic scattering of RayDM against the nucleus 
through two photons exchange. 



the scattering of the WIMP off the external electric field 
generated by the nucleus. In appendix [S] we provide a 
different derivation which leads to a more exact result 
that is valid when the WIMP mass cannot be neglected 
relative to that of the nucleus. At leading order in the 
velocity expansion the amplitude for this process is given 
by 



iM 



iaZ 2 g~ 1 ~ i 




u(k')u(k) u{p')u(p), (24) 



where Z is the nucleus charge, a = 137 -1 is the fine- 
structure constant, Qq is the nuclear coherence scale, 
and the momentum transfer is related to the relative ve- 
locity between the WIMP and the nucleus and the an- 
gle of scattering 9 in the centre of mass frame through 
|q| 2 = 2y?v 2 {\ — cos 8). The function J-{x) decreases 
exponentially for high momentum transfers and is of or- 
der unity near the origin, ^(0) = 2/y / 7r. The spin- 
independent differential cross-section in the centre of 
mass frame is then, 



da 
dcos 9 



2tt 



aZ 2 g 11 Q 




(25) 



Here fi X N is the nucleus- WIMP reduced mass. To a good 
approximation we can use .F(|q 2 | /Qq) ~ -^(0) and so 
the total spin-independent cross-section per nucleon is 
given by 



a 2 Z^g, 
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4tt 2 A 4 A 6 



(26) 



where A is the nucleon number. This is an ex- 
tremely small cross-section for an electroweak scale 
WIMP. For example, taking the nuclear coherence scale 
Qo = ^(O.a + O.SOA 1 / 3 ) -1 , the Raylcigh scale A s = 
500 GeV, and setting g 77 = 1 yields o^ 1 10~ 49 cm 2 
for scattering on xenon. If the Rayleigh scale is consid- 
erably lower, A R < 100 GeV then the scattering rates 



become appreciable. That requires much lighter WIMPs 
than what we set to explore in this work and we leave it 
for a future study to elucidate the prospects associated 
with this part of parameter space. 



IV. COLLIDER PHENOMENOLOGY 

In MiDM, the production mode in colliders is simply 
// — > XX* followed by the prompt decay of the heavier 
state to a photon or a Z boson, as shown in Fig. 7(a)| As 
we emphasized throughout this paper, the phenomenol- 
ogy of MiDM depends crucially on the mass difference 
between the WIMP x an d the heavier state x* an d its 
collider phenomenology is no exception. When the mass 
splitting is small searching for MiDM in collider pro- 
ceeds in the same fashion as searching for other WIMPs, 
namely by looking for mono-jet or mono-photon events 
from an unbalanced initial state radiation (see for ex- 
ample the excellent recent works of [2H426] ). When the 
mass splitting is large the collider signatures of MiDM 
are more similar to that of RayDM. 

In RayDM, since the coupling to WIMPs requires at 
least two vector-bosons, the collider signatures are some- 
what different than usual. There are two different pro- 
cesses that may be searched for: // — » ffxx through a 
vector-boson or photon fusion; and // — > XX^ '■> through 
an intermediate vector-boson, where V = 7, Z, or W ± . 
This last process, shown in Fig. |7(b)| enjoys a larger 
cross-section and results in the production of a photon 
or an electroweak vector-boson in association with large 
missing energy. We therefore concentrate on this pos- 
sibility below. We begin by discussing the MiDM sce- 
nario, followed by RayDM, and finally discuss the actual 
constraints. General formulas for the differential cross- 
sections in the different cases are given in appendix [Cj 



A. MiDM 

In MiDM the dominant mode is the production of 
the heavier state in association with the WIMP through 
// — > 7/Z — > XX* ■ After production, the heavier state 
subsequently decays to the WIMP through the em ission 
of a photon or a vector-boson as shown in Fig. 7(a) 



The differential cross-section for this process is given in 
Eq. ( C-8 ) in the appendix. 



When the splitting between the excited state and the 
ground state is much smaller than the mass AM <C m x 
the resulting photon or vector-boson is too soft to be 
searched for directly. Ref. [57] proposed some interesting 
ideas for looking for iDM in colliders when the mass split- 
ting is in the GeV range and the heavier state decays into 
pions. More generally and without reliance on such spe- 
cialized techniques, the collider phenomenology in this 
case is identical to the usual case of WIMP pair produc- 
tion, but through a dipole operator. It can be searched 
for in a general way by tagging on initial state radia- 



tion. This was nicely worked out in ref. [TS] for mono-jet 
searches in the case of degenerate states (tj x * = m x ) for 
both magnetic as well as electric dipoles. When the split- 
ting is of order the mass and larger, the emitted photon or 
vector-boson may be sufficiently energetic to be searched 
for directly. This can be searched for without reliance 
on initial state radiation hence enjoying a larger cross- 
section. In this case collider searches for mono-photons 
place strong constraints on this scenario as discussed be- 
low in the final part of this section. 

In Fig. [6] we plot the production cross-section for 
MiDM as a function of the WIMP mass for several choices 
of the parameters and the mass splittings. We recall that 
in the case of small splitting where the relic abundance 
is determined by the annihilation into fermions the 77 
signal is independent of the dipole strength whereas the 
collider cross-section scales as the square of the dipole. 
Thus, the cross-section shown in Fig. [6] should be in- 
terpreted as the minimal cross-sections when the dipole 
strength is normalized to yield the correct relic abun- 
dance, /i = ^thermal- For the same reason, the ratio 
o-(pp->7/z->xx ) j s independent of the dipole strength for 

a given mass and choice of /i z / /j, . It is given in Table □ 
and allows for a straightforward comparison between the 
rates in astrophysical processes and the cross-sections in 
colliders. 





Mz/M-, = -tan0 w 


Mz/AS = cote w 


(m x ,AM) 






(130 GeV,0) 


(28 fb) 2 


(0.4 pb) 2 


10" 28 cm^s-l 


lO" 28 cm 3 s-l 


(130 GeV, 100 GeV) 


(16 fb) 2 


(0.2 pb) 2 


lO" 28 cm3 S -! 


10- 27 em's-' 


(300 GeV, 0) 


(1.6 fb) 2 


(23 fb) 2 


lO" 28 em's" 1 


lO" 28 em's- 1 


(300 GeV, 100 GeV) 


(1.1 fb) 2 


(15 fb) 2 


10~ 28 cm 3 s- 1 


10~ 28 em's- 1 



TABLE I. For a given WIMP mass m x and splitting AM = 

m , — m„ the ratio °'( pp ~'" 1 '/ z ~ >x * ) ; s independent of the 
dipole strength and is displayed in this table for several 
choices of the masses and couplings. The production cross- 
section is calculated at leading order for the LHC with = 
7 TeV. 



B. RayDM 

Collider searches for dark matter are inevitably tied 
up in the embedding of dark matter into a complete the- 
ory and RayDM is no exception. In this subsection we 
discuss the possibility of directly producing WIMPs in 
colliders in the RayDM scenario. As we shall see, the phe- 
nomenology is sensitive to the UV physics that resolves 
the Rayleigh operator and as a result is more model de- 
pendent. This is in contrast to the phenomenology asso- 
ciated with direct and indirect detection efforts discussed 
in the previous sections which is insensitive to the details 
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FIG. 6. Minimal production cross-section for MiDM, pp — > 
XX* > through 7/Z in the s-channel against the mass of the 
WIMP. The lower (black) curves show the case of small 
splitting m x « m v „ for two choices of the dipoles' ratio 



Mz/aS 



tan ( 



(solid) and /x z //x = (dotted). The 



dipole strengths in this case are normalized to yield the cor- 
rect relic abundance when the annihilation is dominantly into 
fermion pairs, Eq. Q. While the 77 signal remains invari- 
ant as the dipole strength increases, the collider cross-section 
will increase as fj? . The upper blue curve depicts the case 
when the mass splitting is large m x , — m x = 100 GeV for 
Mz /AS = — tan 9 W . The dipole strength is normalized to yield 
the correct relic abundance when the annihilation is domi- 
nantly into vector-boson pairs, Eqs. ( 9)11 \. 



of the UV physics. It is important to keep this contrast in 
mind when considering the impact of collider constraints 
on the previous sections (this issue also arises in other 
models of DM, see [28 and references therein for some 
related discussion). 

In the case of RayDM, the differential cross-section for 
ff—t XXI at center of mass energy y/s m x through 
an intermediate 7/Z is given by 



1 da 

Ctot d PT 



20 p T 



3/2 



(27) 



where pt is the transverse momentum of the photon in 
the center-of-mass frame, and 



Ctot(// 



XXI) 



aq 2 



3840tt 2 A 6 „ 



(28) 



x (g 2 n + 2g^g lZ v f Z(s) + g 2 z {vj + a 2 f )£ 2 (s)) 



Here qf is the fermion's electric charge, Vf (af) is the 
ratio of its vector (axial-vector) coupling the Z boson to 
its electromagnetic coupling, and £(s) = s/(s — m z ). We 
note that the transverse momentum distribution is such 
that most photons are fairly central, which is important 
for the mono-photon searches. Similar relations can be 
obtained for the production cross-section for // — > xxW . 
In the limit where the WIMP mass and the W-boson 
mass are both much smaller than the incoming center of 
mass energy an identical distribution in px results. This 



motivates mono-W searches, looking for the final state 
W produced in association with the invisible \X pair. 

At LEP for example, where \fs « 200 GeV, one ob- 
tains atot ~ 7 x 10" 3 fb (500 GeV/AJ 6 for 9 X = 0. 
This cross-section is much too low unless the Rayleigh 
scale is brought down considerably. This is in good qual- 
itative agreement with the very thorough investigation of 
ref. [29] where bounds on unparticle production at LEP 
were presented^ 

At Tevatron and LHC, one must convolve the above ex- 
pressions against the parton luminosity functions. The 
resulting cross-section is larger, but for those values of 
the Rayleigh scale where the cross-section is sufficiently 
large to be interesting the theory requires a UV com- 
pletion. One possible UV completion of RayDM is of 
course MiDM in the case when the mass splitting is very 
large. Integrating out the excited state x* one recovers 
the RayDM interactions. So schematically 



MiDM 



-> RayDM. 



(29) 



In this case the Rayleigh scale A R is connected with the 
magnetic dipole \x x through A 3 = m » /2/tiS and one can 
easily translate the results for MiDM in the previous sub- 
section to the case of RayDM. 

Another possible UV completion involves a scalar s and 
a pseudoscalar a, which couple directly to the WIMP as 
well as to the field-strengths of Uy(l) and SU W (2). We 
parametrize this theory with 

C = \d^sd»s - \rn\s 1 + Id^ad^a - \m 2 a a 2 (30) 

cos 6 cos 9 

+ s X X + axisX + -rr^sB^B^ + ^-S-aB^B^ 



4A T 



sin 8 sin 9 



4A T . 



4A,- 



Here we have taken the Yukawa couplings of the scalars 
to the WIMP to be order unity and used A uv to de- 
note the scale of the dimension-5 operators. Integrating 
out the scalars, s and a we generate the different op- 
erators of RayDM. When the scalars are light enough 
to be p roduc ed in colliders the dominant process shown 
in Fig. 7(c) is // — > Vs(a) followed by the decay of s 
(a) to a WIMP pair. We note that this example in fact 
results in a more general version of RayDM where the 
relative coupling of the scalar \X and pseudoscalar X75X 
to the field strengths is arbitrary. This is important in 
the case of comparing direct detection rates (which are 
sensitive to the scalar piece) to indirect detection rates 



8 The Rayleigh operator has scaling dimension of A = 3 in the 
notation of ref. ;29 , which was not considered by the authors 
for good reasons. For A = 2 they find that the unparticle scale 
can be as low as 190 GeV. Attempting to extrapolate to A = 3 
is not very useful since the energy available is greater than the 
cut-off scale and some UV completion is needed to resolve the 
non-renormalizable Rayleigh operator. 
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FIG. 7. The production of WIMPs at colliders is shown in 



Fig. 7(b) for the RayDM scenario. When the Rayleigh scale 
is comparable or lower than the energies involved in t he co lli- 
sion, the Rayleigh operator must be resolved. In Fig. 7(a) we 
show the corresponding process in the case where RayDM is 
the result of integrating out a heavy excited state in MiDM. 
In Fig. ~~ 



Rayleig 



7(c) we show the the process in the case where the 



r operator is resolved in terms of a new scalar. More 



details on the UV completions are provided in the text. 

(which are sensitive to the pseudoscalar piece) . This dis- 
tinction does not play an important role for the purpose 
of collider phenomenology [24- : '2B] . 

In this case, the distribution of the transverse momen- 
tum of the photon in the center of mass frame is given 
by 



1 da 

Otot dp T 



3 p T 



(1- 



1 - 



ip 2 T 



•(31) 



We note that this function is strongly peaked towards 
the kinematical limit pf 1 = ^Q- (l — m 2 /s). This is 
in sharp contrast to typical mono-photon signatures of 
dark matter production where the photon originates from 
initial state radiation and hence its px is dominantly soft. 
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FIG. 8. A plot of the production cross-section for the heavy 
scalars in association with a photon, pp — > ys(a) against the 
annihilation rate of WIMPs to two photons. The mass of 
the scalar is allowed to vary between 2m x and Auv- The 
cut-off scale is set to Auv = TeV, but since both the pro- 
duction cross-section as well as the annihilation rate scale as 
Ayy any other choice can be obtained with a simple rescal- 
ing. The solid (dashed) curves correspond to m x = 130 GeV 
(m x = 300 GeV) whereas the blue (black) curves correspond 
to cos# = 1 (cos# = 0). 



Here the total cross-section is given by 



CTtot (// -> XXl) = 



a 1f 

24A 1JV 



x (g 2 



^1 

s 

2 r„,2 



(32) 
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2g7 7 5 7 z«/C(s) +5 7 z( u / + a f) 6 ( S )) ■ 
The couplings and the function £(s) are defined after 



Eq. (28) above. Similar expressions hold for the axial 
scalar a. In Fig. [8] we plot the production cross-section 
of a photon in association with one of the scalars against 
the annihilation rate of WIMP into two photons. 



C. Limits from Colliders 

Limits on RayDM from colliders come primarily from 
either mono-jet or mono-photon searches. In the case 
of MiDM with a small splitting the production in col- 
liders is observable only through the emission of a gluon 
or photon from the initial state partons. Thus the most 
constraining limits on this scenario come from mono-jet 
searches. The most recent search from CMS [SHj place a 
limit of a few pb in the range m x < 10 3 GeV (see Fig.[9j). 
This is not quite sufficient to exclude the interesting pro- 
duction cross-sections in the MiDM scenario (see Fig.|6|, 
but it comes close. Consulting Tbl.|l]we see for example 
that in the case of hypercharge dominated interactions 
and WIMP mass of m x ~ 100 GeV monojet searches at 
the LHC can begin probing annihilation rates into di- 
photons of the order of 10~ 26 — 10 -27 cm 3 /s. 
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FIG. 9. A plot of the collider constraints on the produc- 
tion cross-section of a WIMP pair coming from CMS mono- 
jet search [30] (solid-red) and CMS mono-photon search [31] 
(dashed-red). The constraints from CDF monojet search |32] 
are only slightly weaker compared with the CMS results. In 
solid black we plot the expected production cross-section in 
the case of MiDM with a small mass splitting and with a 
dipole strength ten times larger than /Uthermal. 



On the other hand, when the splitting is large (AM > 
100 GeV) the photon emitted from the excited state's 
decay is hard enough to be searched for directly. In that 
case the relevant limits are the limits on cr{pp — > p'y) 
coming from mono-photon searches. These are much 
more constraining and are at the level of a(pp — > xXl) ^5 
14 fb for a WIMP mass below the TeV range [3T]. As can 
be seen from Fig. [6] this bound all but excludes MiDM 
with large splittings. In the case of heavy scalar produc- 
tion in association with a photon one can again compare 
the production rate directly to the limits on a(pp — > $7) 
since the resulting photon is hard. Since the production 
cross-section in the case of scalars is generically lower 
(see Fig. [8]) the current constraints are not quite strong 
enough to exclude this scenario, but they are now probing 
the most interesting parts of parameter space. 



V. CONCLUSIONS 

The effective theory describing the interactions of a 
Majorana WIMP with photons is of critical importance, 
given that our best indirect detection searches come 
through monoenergetic 7-ray lines, and direct detection 
is clearly sensitive to scattering through a photon ex- 
change. Interestingly, this effective theory is quite re- 
stricted: in the presence of a nearby excited state, there 
is the possibility of an interaction with electromagnetism 
via a dipole transition to the excited state (or Magnetic 
Inelastic Dark Matter or MiDM); in the absence of a 
nearby state, the leading operator comes in the form 
XXW^W^ or xxBfivB^" or its the pseudoscalar and 
CP violating equivalents. These two scenarios have re- 
lated, but distinct phenomenology. 



Remarkably, in the case of MiDM both the size of the 
signal in direct detection and 77+7Z signatures are in- 
dependent of the size of the dipole, with the relic abun- 
dance suppression precisely canceling out against the en- 
hanced scattering and annihilation cross sections. This 
offers a surprising concordance whereby the annihilation 
rates into 77 is in the range to explain the tentative ex- 
cess in gamma rays at around 130 GeV and possibly ex- 
plain the DAMA annual modulation. MiDM predicts a 
secondary line at around 114 GeV from 7Z with a rel- 
ative rate of about 1 : 3 compared with the 77 line at 
130. Fermi should be able to test the 7 ray signature 
and, for small mass splitting m , — m x w 100 keV, the 
MiDM scenario also predicts collision rates with nuclei 
that can now be tested at direct detection experiments. 
The production rates in colliders are below the current 
sensitivity of the LHC for thermal cross sections, but can 
exclude some regions of parameter space where this par- 
ticle constitutes only a fraction of the total dark matter. 
The case of the thermal WIMP constituting all of the 
dark matter should be observable in the near future. We 
showed that given the scaling of the different quantities 
involved, the concordance is in fact independent of the 
dipole strength and is maintained even with an increased 
dipole strength where MiDM forms only a fraction of the 
total DM. Model-dependent corrections outside of the ef- 
fective theory can change this result, however. Moreover, 
constraints from colliders place an ultimate limit on such 
an increase in the dipole to be no more than O(10). 

In the case of RayDM it is also possible to simultane- 
ously achieve the right relic abundance as well as rates in 
the range now explored by gamma ray observations. But, 
in contrast with MiDM, it favors stronger coupling to the 
SU W (2) vector-bosons than to hypercharge. If RayDm is 
describing only the interactions with photons, however, 
and freezes out through some other channel, coupling to 
hypercharge alone gives a good description of the data. 
RayDM predicts a ratio of 7Z to 77 of 1 : 5 when cou- 
pling to hypercharge dominates, or about 5 : 2 in the 
more likely case of dominant coupling to SU W (2). Un- 
fortunately, the direct detection prospects in this case 
are gloomy as the collision rates with nuclei due to two 
photon exchange are much too small. In contrast, this 
scenario offers interesting phenomenology in colliders in- 
cluding mono-photon, mono-Z, and mono-W signatures 
with rates that can now be probed at the LHC. 

There are a few interesting variations on the scenar- 
ios we have discussed. A particularly natural scenario is 
MiDM+RayDM, where the 77 signal is naturally boosted 
in an MiDM model by the presence of an additional hy- 
percharge Rayleigh operator. Such an operator is gen- 
erally present and would be expected to often dominate 
the 77 signal from these models. 

An alternative possibility is that some amount of 
hypercharge-dominated RayDM is just a subdominant 
component of the dark matter. Since the density scales 
as p ~ (o"u)~ nn , the overall rate scales as p 2 {w) ann ~ 
( av )ann- Thus, rather than having all dark matter 
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annihilate to 77 with a cross section (<Jv) ann ~ 3 x 
10~ 27 cm 3 s _1 , we could have a cross section ~ 10 x 
( av ) thermal an d yield the claimed 77 signal from a sub- 
dominant component of dark matter. 

While a number of opportunities exist to distinguish 
the MiDM scenario from a RayDM scenario, there is 
another important difference: in RayDM, in particular 
when the dominant operator is XX W^vW^", there is a 
sizable hadronic annihilation channel (via W's and Z's) 
compared to 77. In contrast, for MiDM, the // channel 
is not present in the late universe as it is only present for 
X*X annihilations rather than XX- Limits on the contin- 
uum photon emissions such as those from dwarf galaxies 
[33] [34] or the galactic center [35H37] could potentially 
distinguish these scenarios. 

Ultimately, if a Majorana dark matter interacts signif- 
icantly with light, there are a number of conclusions that 



can be drawn right away. While direct detection signals 
require a nearby state, collider signatures do not. The 
era of data - approaching dark matter with direct, in- 
direct and collider experiments, may be on the verge of 
revealing its nature. 
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Appendix A: Elastic Scattering due to two photons exchange 

In this appendix we derive the amplitude for the elastic scattering of the WIMP on the nucleus due to two photons 
exchange, Eq. (24 1. We begin by noting that there are several separate scales in the problem: A R , the high scale 
associated with the Rayleigh operator; m N , and m , the masses of the nucleus and the WIMP, respectively; q 2 , 
the momentum exchange, approximately 100 MeV; R N , the nuclear coherence size approximately 100 MeV; q°, En, 
fiv 2 - the kinetic energies involved, approximately 10 keV. By itself this diagram is logarithmically divergent, but 
inclusion of the charge form- factor provides a natural cut-off at a momentum scale around R~ . Therefore, the 
momenta running in the loop are non-relativistic and we evaluate this diagram using known techniques from heavy 
quark effective theory [38] - In appendix [B] we show that using this technique one can recover the results obtained in 
ref . [35] where the second order Born cross-section was used to calculate the elastic channel of usual iDM [?D] . We 
begin by writing the momentum of the intermediate nucleus in the usual velocity expansion, 

P = m N v + P, (A-l) 

where v = (1, v) is the 4- velocity of the nucleus. The propagator for the nucleus can be approximated as, 

f + m N 1 l + ji 



p2 



(A-2) 



% P-v + i0 2 

where (1 + f)/2 is the projector onto the two large components of the 4-spinoir] The projector causes the QED vertex 
of the nucleus to simplify to iZev^ instead of the usual iZe^. This embodies the fact that in the non-relativistic 
limit, the polarization of charged particles does not change under the exchange of a photon. Given the above, the 
amplitude associated with this diagram at leading order in the velocity is given by 



iM = (u(k') u{k)) v( 9 2 ) (u(p') «(?)), 



with 



and 



9jf 



dH 



i + i> 



-9m 



(A-4) 



(A-5) 



-9u 



A3 J (2tt)4 \{l + q /2) 2 ) VH + 9/2) 2 
x ({-I 2 + g 2 /4) -{1 + q/2)^ {-I + q/2f' 



(P + 1 + q/2) - v + i0' 



(A-6) 



9 In the diagram above, P 
can write, 



p + / + q/2 and since P 



q we 



1 



1 



(P- 



-<?/2) 2 



(p + I + q/2) -v + i0 2 



(A-3) 
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Here I^^i^ 2 ) nas dimensions of inverse mass square. Since we are interested only in the leading order in the velocity 
only the time-like indices are important and we obtain 

e 2 Z 2 g, n f dH l 2 -q 2 /4 



A 3 



x F 



F( 



1 



(A-7) 



(2tt) 4 D 

Here, bold face letters denote 3-vectors, and we included the charge form-factor F(|q|) by hand. The denominator is 

D = (l + q/2) 2 (-l + q/2) 2 (p° + l Q + q°/2 + iO). (A-8) 

In order to allow for exact evaluation of this integral we choose to work with the Helm form-factor, which is a function 
of the momentum exchange q 



F(|q|) =e- q2 / Q o. 

Here Q = \Z6R~ 1 , and the nuclear radius is (see for example the excellent review by Salati [4"I] ) 

R„ = fm X ( 0.3 + 0.89A 1/3 ^ 



Neglecting the overall constant in front, the integral of Eq. (A-7) can now be written as 

dH l 2 -q 2 /4 / l 2 + q 2 /4 s 



ho (v ) <* 



(2tt)' 



D 



x exp 



Ql 



(A-9) 
(A-10) 

(A-ll) 



using the usual Feynman parameter together with an HQET parameter with the dimensions of energy 



The denominator D, given in Eq. (A-8), is a factor of several separate propagators that can be combined together 

(A-12) 



1 

D 



d:r 



d£- 



{£ (p° + 1° + q°/2) + l 2 + q 2 /4 + (1 - 2x)l ■ q + lef ' 
Shifting the momentum variable — > + h ((1 — 2x)q^ + Sgo^) we can write the integral as 



dxd£- 



d-n 



dl° 



N 



(2tt) 3 J 2tt ((;o)2 _ A + ie y 



exp 



l 2 - (1 - 2x)\ ■ q + (2 - 4x + 4x 2 )q 2 /4 



Ql 



where 



N = 2 (l 2 - (1 - 2x)\ ■ q + x(x - l)q 2 ) , 
A = l 2 + ^ + x(x - l)q 2 -£(p° 2 + xq°) 



(A-13) 
(A-14) 



The mixing term in A can be neglected as it is always much smaller than the other two terms. Either way, the integral 
over Z° can be done exactly and yields 



00 dl° 



3/ 



-oo 2 ^ ((l°) 2 - A + ie)- 



16A 5 /2- 



(A-15) 



The dependence of the form-factor on 1 • q = |1| |q| cos 9 causes the integration over d 3 l to be slightly more complicated 
than usual. We proceed by first doing the integral over cos 9 followed by the integral over the dimensional Feynman 
parameter £. We are left with two integrals, one over the Feynman parameter x and the other over the radial 
component of the spatial momentum |1|. Defining the dimensionless variables I = |1| /Qq and q = |q| /Qq we arrive at 
the result quoted in the text in Eq. ( 24 ) 



iM 



iaZ 2 g 11 Qq 
4 AT 



T 



u(k')u(k) 



where, 



^(9) = ~ 

7T 



dx 



dl ^r- 



P 



(I 2 + (1 - x)xq 2 ) 2 
x cosh ( (1 — 2x)lq 



x exp (^-l 2 - q 2 (\ - x + x 2 ) 
I 2 - (1 - x)xq 2 + 1 



(A-16) 



(A-17) 



(1 - 2x)lq 



sin I (1 — 2x)lq 
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Appendix B: Derivation of the elastic channel of iDM using HQET 

The inelastic Dark Matter scenario of Ref. [40] involves a dark matter state x thai interacts with normal matter 
only through a transition involving an excited state x* separated in mass by AM. Ref. [39] considered the particular 
case where the interaction with the SM is through a new massive U'(l) vector-boson that kinetically mixes with 
hypercharge, a Holdom boson [35]. They computed for the first time the contribution to the elastic channel from the 



second-order diagram shown in Fig. 10 in the case where the mass splitting is much larger than the kinetic energy 
available, AM > MeV. In what follows, we reproduce this result using the HQET methods discussed in the text and 
appendix [X] For simplicity we ignore the charge form-factor. 

Writing the four-momentum of the WIMP as k = m x u + k with u = (1, u) being the 4- velocity, we can express the 
momentum of the excited state as, 

K = k-l-q/2 = m x u + k-l-q/2. (B-l) 

Since the mass splitting, AM is greater than the kinetic energy, the fermionic propagator of the excited state can be 
written as, 

* + m x - „ 1 i±if mo) 

K 2 -m?, ~ AM 2 { ' 



The propagator for the nucleus is as given in Eq. ( A-2 ) . Replacing 7^ by four- velocities u M in the vertices and 



considering the leading order term in the velocity expansion the amplitude is given by 

iM = (u(k') u(k)) I(q 2 ) (u tt(p)) , (B-3) 

with, 

r, „ 167TVZW f dH 1 

J < g2 > = 2x AM J(2^D- 

Here the factor of 2 arises from a second diagram where the A' lines cross and which contributes equally at this 
order. Here K is the kinetic mixing parameter, a' — e 2 /4tt is the U'(l) charge. D is the same denominator as 



previously considered and presented in Eq. ( A-8 ) except that the massless photon propagators are replaced by a 
massive propagator (q 2 + ie)^ 1 — > (q 2 — m 2 A , + ie)^ 1 . Following the same steps as in appendix \M we introduce the 
Feynman parameters x and £ to combine the denominator, shift the loop momentum Z M — > ^ J\l — 2x)q ll + Sgo^), 
and then integrate over dl°, 

i\ AM J J, J J (2i)'(12 + A /)V2 

Here A' = £ 2 /4 + x(x — l)q 2 + m 2 A , — £(p° + xq°). Neglecting the term linear in £ as before, all the integrals can 
easily be done to yield 

i ( \&ir 2 K 2 Z 2 aa'\ Arctan (|q|/2mA') m A ,^o 2i: 2 Z 2 aa / K 2 
m) = -^{ AM J H |q|AM ■ (B " 6) 

The differential cross-section in the CM frame is simply 

d(J U 2 r n„ , s 4a 2 a'2Z 4 K 4 u 2 , / Iql \ 2 m A ,-yO Tr 2 aa'Z 2 K 2 u 2 

— = i-^I 2 ( q ) = „ . - r Arctan „ A w , B-7) 

4tt 2 q 2 AM V 2m A'/ q 2 AM ' v ; 

which reproduces the result obtained in ref. [33] in the appropriate limit. 



Appendix C: Differential Cross-Section Formulas for MiDM and RayDM Production at Colliders 

In this appendix we give the formulas for the production cross-section and distributions for the process // — > xxY 
where V — ^f,Z,W ± through the Rayleigh operators of Eq. (3| as well as related processes in UV completions of 
RayDM. A convenient way of presenting the results is obtained by treating the XX system as having four-momentum 
Pi and mass p 2 . The 2 — > 3 phase-space factor can then be written as 

dpi 

dPS 3 (pfpj ^Pv+Px+ Px) = dPS 2 (PfPf -> Pv + Pi) x — ^ x dPS 2 (Pi -^Px+ Px) ■ (C- 1 ) 
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FIG. 10. The second Born amplitude for the process \N — > \N through an intermediate excited state x* ■ A second diagram 
with the A' lines crossed contributes the same as this one. 



Averaging (summing) over initial (final) polarization for the process shown in Fig. |7(b)[ and neglecting the incoming 
particles' masses, the matrix element squared is 



ly w ,2 = (gg +9 2 a)9 2 w spj 1 



(C-2) 



pol 



o| (l - W/pj) (m 4 v + (s - pi) 2 + 2m 2 v (s - t) + 2(s - p 2 )t + 2t 2 ) 
+a 2 A [m% + {s- p 2 4 ) 2 + 2(s - p\)t + 2t 2 - 2m 2 v {s + <)) 



2f 

where g v>a are the vector and axial couplings of V to //, gvv arc defined in Eq. and the Mandelstam variables 

s = (pf +Pj) , t — (pf — pv) 2 , and u — {pf — Pa) 2 are defined as usual [33]. Here we kept explicit the separate 

contributions from the scalar XX^F (axial XlsXFF) piece by preceding it with as (cia)- This separation will prove 
useful below when discussing the corresponding formulas in the UV completion of RayDM with heavy scalars. The 
integral over dPS2 (p4 — > p^ + p x ) is straightforward and can be done in its entirety since the squared amplitude in 
Eq. (C-2 1 contains no dependence on the XX system's angular distribution. The integral over the azimuthal angle of 
dPS2 {PfPf Pv +P4) can also be done and the differential cross-section is then given by 



d 2 a(ff^xxV) 



{9l 



9 2 a)9 2 V v Pi 



dcosd dp: 



2048tt 3 A6 ( s - m 2 v f 




(l - Am 2 Jpl) (ru* + (p 2 - s) 2 + 2m 2 v ((^^) s -pi)) 



+a 2 A (my + (s - Plf 



2m 2 v (p 2 4 + s)) 



(C-3) 



where 9 is the angle between the incoming fermion / and the vector V in the center-of-mass frame, and X(x, y, z) 



y 



z — 2xy — 2yz — 2zx is the usual kinematic function. The angular region is — 1 < cos 9 < 1 and the mass 



of the XX system is in the range 4m 2 < p\ < (y/s 
d 2 v(ff^XXl) _ag^ L p 2 (l 

5127T 2 



d cos 9 dp 2 



cos 



A 6 



The photon case is particularly simple and yields 

3/2 



4m; 



(a 2 A + a 2 (l-4m 2 x /p 2 )) 



1 



Pi 



(C-4) 



This can be integrated exactly to yield the total cross-section that is quoted in Eq. ( 28 1 for the case where y/s ^> m 



What is often of more interest is the transverse momentum distribution of the vector boson. The transverse momentum 
in the centre-of-mass frame is given by px = Pv sin 9 and the momentum of the vector in the center-of-mass frame is 

pv = ^pA 1 / 2 ^1, ^j 1 , . This can be used to obtain the differential distribution 



da 
dp T 



dpi dcos9 



d cos 9dp 



2 5 (p T (cos 9, pi) -p T ) 



(C-5) 
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Needless to say, for hadronic colliders such as the LHC and the Tevatron these expressions have to be convolved 
against the appropriate parton distribution functions. We have verified the formulas above against the Madgraph 4 
package [4"4"j . 

Similarly, in the case of heavy scalars, by putting the \X momentum on-shcll p\ = m 2 we obtain the matrix-element 
squared for the process // — > js through 7/Z 

zEl^| 2 = ! S^ t - AjL * (g^ + 2g, n g, z v f a S ) + + a/) a £ 2 (*)). (C-6) 

pol u v 

Here s,t, and u are the usual Mandelstam variables. Use of the relation 

da [ Pt 1 \ do 



dpr \ P yjp 1 - p\ J d cos 6> ' 



(C-7) 



leads to the differential cross-section given by Eq. (31 1 . Here P — 2 <rj ^ s ^ ne momentum of the photon in the 

center of mass frame. Similar expressions hold for the other distributions involving Z and W . 



Finally, in the case of MiDM, the production cross-section is given by 

/' ' ' ' J ' /' 



^ 1 + 2vAti{s) + (vj + a}Ae(s) ] (C-8) 



d cos 6 



\ 



( m 2 m 2 \ ( AM 2 \ ( 9 (m +m„) 2 . 
A I 1, J x (1 ) [ sin 2 9 + ^- (1 + cos 2 ■ 



where 9 is the scattering angle in the centre of mass frame, and w 7i z ( a 7,z ) is the vector (axial) coupling of the 
corresponding vector-boson to the incoming fermions. 
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